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Abstract 

Renormalization group methods can be applied to the nuclear many-body problem 
using the approach proposed by Shankar. We start with the two-body low mo- 
mentum interaction V^ow k and use the RG flow from the particle-hole channels to 
calculate the full scattering amplitude in the vicinity of the Fermi surface. This 
is a new straightforward approach to the many-body problem which is applicable 
also to condensed matter systems without long-range interactions, such as liquid 
^He. We derive the one-loop renormalization group equations for the quasiparticle 
interaction and the scattering amplitude at zero temperature. The RG presents an 
elegant method to maintain all momentum scales and preserve the antisymmetry 
of the scattering amplitude. As a first application we solve the RG equations for 
neutron matter. The resulting quasiparticle interaction includes effects due to the 
polarization of the medium, the so-called induced interaction of Babu and Brown. 
We present results for the Fermi liquid parameters, the equation of state of neutron 
matter and the ^Sq superfluid pairing gap. 
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1 Introduction 



Fermi liquid theory is a prototype effective theory, invented by Landau in the 
late 1950's [1,2,3]. In this theory the properties of normal Fermi liquids at zero 
temperature are encoded in a few parameters, which are related to the effec- 
tive two-fermion interaction. Over the years, Fermi liquid theory has proven 
to be an extremely useful tool for studying normal Fermi liquids, in partic- 
ular hquid ^He [4]. The main assumption is that the elementary, low-lying 
excitations of the interacting system are relatively long-lived quasiparticles, 
with properties resembling those of free particles. The remaining strength in 
the spectral function is distributed over modes that add up incoherently. In 
the language of the renormalization group, the incoherent background is in- 
tegrated out into the quasiparticle interaction, which can be determined by 
comparison with experiment or by microscopic calculations. 

Babu and Brown later realized [5] that, in order to satisfy the Pauli principle 
in microscopic calculations, one has to take into account not only the particle- 
hole channel considered by Landau, which gives rise to the propagation of zero 
sound, but also the exchange channel thereof. This is taken into account in 
the induced interaction, which incorporates the response of the many-body 
medium to the presence of the quasiparticles. 

Almost a decade ago, Shankar revived the interest in Fermi liquid theory by de- 
veloping a renormalization group (RG) approach to interacting Fermi systems 
(for an introduction see Shankar [6] and the lecture notes of Polchinski [7]). 
Since Fermi liquid theory is restricted to low-lying excitations in the vicinity of 
the Fermi surface, normal Fermi systems are amenable to the renormalization 
group, where one has precisely such a separation of modes in mind. More- 
over, by taking the loop contributions from both particle-hole channels into 
account, the RG flow remains antisymmetric, i.e., at any step of the renormal- 
ization, the scattering amplitude obeys the Pauli principle and consequently 
the solution obeys the Pauli principle sum rules. Thus, the antisymmetry of 
the scattering amplitude, which was the guiding principle of Babu and Brown, 
is realized naturally in the RG approach. Recently, we obtained additional RG 
invariant constraints [8] . These approximate constraints are analogous to those 
obtained by Bedell and Ainsworth for paramagnetic Fermi liquids [9]. 

In this paper, we solve the one-loop RG equations in the particle-hole channels 
at zero temperature and in three dimensions. We work in the approximation 
that both particle-hole momentum transfers are small compared to the Fermi 
momentum. The justification for this approximation is two-fold: In Fermi liq- 
uid theory the long-wavelength excitations play a central role. Hence, our 
primary aim is to treat these correctly. Furthermore, for the case of nuclear or 
neutron matter, the dependence of the induced interaction on Landau angle is 
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fairly weak [10,11]. Thus, we expect an expansion in momentum transfers to be 
accurate at least for the Fermi liquid parameters.^ In this approximation the 
treatment of the particle-hole phase space with cutoffs is enormously simpli- 
fied. An RG analysis of a schematic model, including the complete particle-hole 
phase space in two dimensions, is discussed in [12]. Improved RG equations, 
beyond one- loop order, are presented in [8]. 

Here we apply the RG techniques to neutron matter, where the tensor force 
does not contribute in the S-wave. In this system, a further physical motivation 
for expanding in small momentum transfers over the Fermi momentum is pro- 
vided by the unique low momentum nucleon-nucleon interaction Viow k [13,14], 
which we employ as the starting point for the RG flow. As shown in Figs. 2 
and 4 of Bogner et al. [14], ^owk is significant only for relative momenta 
k < 1.3 fm~^. Moreover, from the evolution of V^owk shown in Fig. 5 of [14] 
we deduce that the important momentum modes are around the pion mass 
m^r ~ 0.7 fm~^, since higher momenta do not renormalize V^owk considerably. 
Therefore, the typical momentum transfer is small compared to the Fermi 
momentum of neutron matter at nuclear matter density kp = 1.7fm~^. 

The main idea of the RG approach to the many-body problem is to "adiabati- 
cally" include the in medium corrections to the effective interaction by solving 
the RG flow equations in the relevant channels. In this exploratory calculation 
we include the particle-hole channels, which play a special role in Fermi liq- 
uid theory. The main effect of scattering in the particle-particle channel, the 
removal of the short-range repulsion, is taken care of by using V^ow k as the 
starting point for the RG flow. The low-lying excitations in this channel, which 
are responsible e.g., for superfluidity, are not included. These are then treated 
explicitly in Section 4.2, where we compute the superfluid gap by employing 
BCS theory for the particle-hole reducible scattering amplitude. 

In the RG approach, the width of the momentum shell that is integrated out 
in one iteration is a "small parameter" . The full scattering amplitude is an RG 
invariant quantity. In a conventional approach, this is calculated from a "set 
of diagrams" while in the RG approach it is continuously evolved from the 
starting interaction by gradually including the many-body corrections from 
the narrow momentum shells one at a time. In this sense the RG provides 
a method for dealing with strongly interacting systems, where perturbation 
theory fails. We stress that this application of the renormalization group to 
many-body systems is only indirectly related to the well known RG treatment 
of critical phenomena. 

We derive the one-loop renormalization group equations for the quasiparticle 
interaction and the scattering amplitude at zero temperature. The evolution of 

This may change when effects of the tensor interaction are included in the induced 
interaction. 
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the effective mass is included in the RG flow, as well as a simplified treatment 
of the renormalization of the quasiparticle strength. In the forward scattering 
limit, the role of the energy transfer is in the RG taken over by the cutoff in 
momentum space. We find that, in the long-wavelength limit, the dependence 
of the effective four-point vertex on A/g corresponds to the behavior with uj/q 
in the microscopic derivation of Fermi liquid theory by Landau. 

Finally, we present the solution of the RG equations for neutron matter. These 
are obtained by employing the unique low momentum nucleon-nucleon inter- 
action as the initial condition of the RG. Our results include the Fermi liquid 
parameters in the density range of interest for neutron stars as well as the 
full scattering amplitude for general (non-forward) scattering on the Fermi 
surface. Using the Fermi liquid parameters, we compute the equation of state, 
including polarization effects, by integrating the incompressibility [15]. Finally, 
we compute the ^Sq superfluid pairing gap using weak coupling BCS theory. 
This is an application of our approach, which probes the angular dependence 
of the scattering amplitude. We generally find very good agreement with the 
results obtained in the polarization potential model by Ainsworth, Wambach 
and Pines [16,17]. However, it is worth noting that much of the model depen- 
dence, inherent in the work of Ainsworth et al, can be avoided in the RG 
approach. 



2 Renormalization group at one-loop 

Shankar suggested that the relevant modes of a normal Fermi system may 
be isolated in a similar manner as for critical phenomena. In analogy with 
a Wilson-Kadanoff treatment of the latter [18], he proposed to separate the 
"slow" modes from the "fast"[^ ones by imposing a cutoff around the Fermi 
surface [19]. Using a loop expansion, one can systematically generate the ef- 
fective quasiparticle scattering amplitude and the effective quasiparticle in- 
teraction among slow quasiparticles and slow quasiholes at a given scale A.[^ 
It has been shown that, in this approach, Fermi liquid theory emerges in the 
infrared limit when the cutoff is taken to zero [6,7,19,20,21,22,23,24,25]. 

Recently we derived the flow equations for the effective scattering amplitude 
and the quasiparticle interaction from the induced interaction of Babu and 
Brown [8]. In this paper, we employ the one-loop RG equations to compute 

^ With "slow" we refer to the modes within a shell of width 2A centered around 
the Fermi momentum kp, whereas the "fast" modes lie outside this shell. 
® The effective scattering ampHtude at the scale A interpolates between the two- 
body irreducible interaction and the full scattering amplitude as described in the 
introduction. 
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the full scattering amplitude and the quasiparticle interaction on the Fermi 
surface starting from the vacuum two-body interaction. As we decimate down 
to the Fermi surface by letting A — 0, we obtain not only the forward scat- 
tering amplitude for low-lying quasiparticle-quasihole excitations, but also the 
amplitude for general (non-forward) scattering processes of quasiparticles on 
the Fermi surface in the small momentum approximation. 

In a loop expansion, the induced interaction corresponds to a summation of 
the direct and exchange particle-hole channels in a particular kinematic win- 
dow, while loop contributions from the particle-particle (BCS) channel, are not 
summed by the RG equations derived in [8]. Thus, the contributions from the 
BCS channel should be included in the particle-hole irreducible driving term. 
In practice, one usually restricts oneself to the two-particle ladders with the 
bare nucleon-nucleon interaction (the Brueckner G matrix) or a simple gener- 
alization thereof, where the lowest order renormalization of the quasiparticle 
strength is included [10]. 

In the present work, we explore the one-loop truncation of the RG equations. 
When only one channel is considered, the one-loop RG equation is exact in 
the sense that it is equivalent to the corresponding scattering equation. When 
both the direct and exchange channels are included, the scattering amplitude 
remains antisymmetric under the RG flow. Calculations with the induced in- 
teraction for liquid ^He and nuclear matter show that the systematic inclusion 
of exchange channel contributions is crucial for obtaining a realistic descrip- 
tion of the effective interaction [9,10,11,26,27,28,29]. Therefore, we expect the 
one-loop RG to capture the essential physics of the soft modes near the Fermi 
surface. In this approximation, our RG equations are equivalent to those ob- 
tained previously by Dupuis [20]. 

The two one-loop particle-hole contributions to the effective four-point func- 
tion are shown in Fig. 1, where the integration is restricted to a fast particle 
and a fast hole. The effective four-point vertex 7 depends on three momentum 
variables, the cutoff and the energy transfer u in the zero sound channel,^ 

7 = 7((^,q),q',p + p';A), (1) 

where q' = p — p'- Under exchange, the momentum transfer q and q' have 
to be interchanged, whereas the particle-pair momentum p + p' remains in- 
variant. We note that in the limit |q| 0, the scattering amplitude describes 
long- wavelength excitations in the zero sound channel, while |p + p'| — ^ 

''' For energy- independent initial conditions, one can set the energy transfers lo 
and Lo' to zero in the RG equations. Here we retain only the lo for the purpose 
of discussing the limits uj/q ^ and q/uj — > 0. The role of uj for the purpose of 
obtaining the forward scattering amplitude and the quasiparticle interaction is then 
taken over by the cutoff A, as we will show in the next section. 
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Fig. 1. The one-loop contributions to the RG equation. (A): The zero sound (ZS) 
channel with momentum transfer q and (B): the exchange channel (ZS') with mo- 
mentum transfer q'. (Our convention for the momentum flow is in the direction 
of the particle and hole arrows.) As defined in the text, 7 denotes the effective 
four-point vertex. 



corresponds to the BCS pairing singularity. 

We define the effective quasiparticle interaction /(q'; A) at the scale A by the 
limit 



and the effective quasiparticle scattering amplitude, extended to finite mo- 
mentum transfer q, by 



In order to allow a proper treatment of the momentum dependence in the 
exchange channel, we generalize the effective interaction /(q'; A) to finite q, 
/(q, q';A), keeping in mind that the Landau quasiparticle interaction is ob- 
tained for |q| =0. 



/(q'; A) = lim 7((u;, q = 0), q', p + p'; A), 

OJ — >\J 



(2) 



a(q, q'; A) = 7((^ = 0, q), q', p + p'; A). 



(3) 
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The RG equations at one-loop and zero temperature are given by [8] 



fast.A 
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— /(q,q';A) = -^. — g 
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/ / P + P' , q // *\ / / // P + f 



(4) 



(5) 



where g denotes the spin-isospin degeneracy, for neutron matter g = 2, and 
is the quasiparticle strength at the Fermi surface. The effective scattering 
amphtude includes contributions from both the ZS and ZS' channel, whereas 
for the effective interaction only the exchange contribution, which corresponds 
to the induced interaction [5], remains. The vertices in the one-loop diagrams 
are given by the scattering amplitude at the current scale. Consequently, the 
renormalization due to the fast modes in both particle-hole channels is in- 
cluded. Furthermore, since we will employ the unique, energy-independent 
low momentum nucleon-nucleon interaction of Bogner et al. [13,14] as initial 
condition of the RG flow, we have set the energy transfer to zero, in particular 
a; = 0, in the above RG equations. Note that the derivatives with respect to A 
on the right hand side of Eqs. (4,5) act only on the phase space factors. Fur- 
thermore, in order not to overload the equations, we have suppressed the spin 
variables. They will be discussed in detail below. It is apparent that the RG 
equation, Eq. (4), preserves the antisymmetry of the scattering amplitude un- 
der the RG flow. Hence, the Pauli principle sum rules [3,30] are automatically 
fulfilled for the resulting Fermi liquid parameters. 



As motivated in the introduction, we evaluate the one-loop beta functions in 
the long-wavelength approximation. Obviously, the phase space for the inter- 
mediate particle-hole pair in the ZS channel is open for A < g/2, where q 
denotes |q|, with an analogous expression for the ZS' channel, see Fig. 2. For 
small g, q\ the integration over p" in the beta functions reduces to an integra- 
tion over a solid angle on the Fermi surface. We find for the RG equations at 
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Fig. 2. The phase-space factors in the RG equation, Eq. (4), are non-zero in the 
fihed areas. The dashed circles denote the cutoff around the Fermi surfaces. 

zero temperature^ 

A ^a(q, q'; A) = Q{q - 2A) Pzs[a, q, A] - Q{q' - 2 A) /5zs'[«, q', A] (6) 

2A 

= Q{q - 2A) zl — N{A) 



rffip// sign (cos (Vc 



2A/g<|cos(ep„q)|<l An COs(6'p//q) 



/ p + p , q " \\ f // p + p , q A 
X a(q, — — + - p ; Aj a(q, p — + 

- e{q' - 2A) zl — {q ^ q'} (7) 



A ^/(q, q'; A) = -Q{q' - 2A) zl ^ N{A) 



X 



dVlpn sign(cos(Vc 



2A/g'<lcos(ep„q,)l<l 47r COs(^^p"q' 

/P + p' , q „//. a\ „// p + p' , q. 

2 



X a(q,— — + --P ;Aj a(q,p — + -;A 



The factor A^(A) = g m*(A) fcp/ 2 vr^ denotes the density of states at the Fermi 
surface. The effective mass is determined self-consistently from the effective 



^ Contrary to statements in the hterature, the RG equations at zero temperature 
are very transparent. The 6 functions appearing in the beta function can be handled 
in a straightforward manner and the phase space integrals can in fact be evaluated 
analytically. A non-zero temperature introduces a further scale, which competes 
with the cutoff and consequently must be treated with special care. 
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interaction 

m*(A) 1 



m ^ ornkp 



(9) 



67r2 

This enforces Gahlean invariance on the system. The effective mass is continu- 
ously adapted as the flow equations are integrated towards the Fermi surface. 
In Eq. (9), /i(A) is given by the projection of the quasiparticle interaction 
on the / = 1 Legendre polynomial. If one restricts the analysis to the ZS 
channel and assumes that the Fermi liquid parameters are given, the solu- 
tion of the flow equation reproduces the relations between the quasiparticle 
interaction and the forward scattering amplitude. In the appendix we give a 
detailed derivation of the beta function in Eqs. (6-8) and the solution of the 
RG equation in the ZS channel. 

For A ^ k-p, the particle-hole contributions vanish in the momentum range 
of interest q,q' ^ 2kp (see the explicit expressions given in the appendix). 
Hence, we can start the decimation at Aq = k-p with a particle-hole irreducible 
boundary condition on the RG equations, Eqs. (7,8). This boundary condition 
may be identified with the direct interaction of [5] . We approximate the direct 
interaction with the vacuum two-body low momentum interaction Viow k at a 
cutoff scale Ayj^^ = ^/2kp, which will be motivated in Section 4. This choice of 
the cutoff approximately accounts for the Pauli blocking of intermediate states 
in the BCS channel, which suppresses large structures, like the quasi bound 
state in the ^Sq channel and tames the short-range repulsion of the nucleon- 
nucleon interaction. Since, in isotriplet channels, Vio^, k is almost independent 
of the cutoff over a large range, the precise value of A^j^^ ^ is not crucial. 
Furthermore, we allow for a minimal in medium correction of the starting 
values, in terms of the renormalization of the quasiparticle strength Zkp- To a 
certain extent this accounts for the contribution of higher order particle-hole 
irreducible diagrams to the direct interaction [10]. At the initial cutoff Aq, we 
thus have 

a(q, q'; Aq) = /(q, q'; Aq) = 7vacuum(q, q'), (10) 
where 7vacuum(q; q') will be specified below. We solve the RG fiow with two 
different assumptions for the Zkp factor. In one case, we will use a static, 
density-independent mean value of zl^ = 0.9, which remains unchanged under 
the RG. In the other case, we compute the Zkp factor dynamically, by assuming 
that the change of the effective mass from the initial one, owing to the direct 
interaction, is due to the Zkp factor alone. We then adjust the factor z^^ in 
the RG equations, Eqs. (7,8), self-consistently as we decimate to the Fermi 
surface, i.e., we start with Zkp{AQ) ~ 1 and change the z^p factor according to 

1 _ ^*(A) ^^^^ 



Zkp{A) m*{Ao)' 

In this approximation we account for the momentum dependence of the self- 
energy corresponding to the direct interaction in m*(Ao) and make the rea- 
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p'+q 




Fig. 3. The configuration of quasiparticles and quasilioles on the Fermi surface. The 
arcs denote the angles 9^ and 6^' opposite of q and q', respectively. 

sonable assumption that the modification of the effective mass due to the 
induced interaction is dominated by the energy dependence of the self en- 
ergy, the so called E-mass [31]. We shall see later that this approximation 
leads to a Zkp factor in rather good agreement with Brueckner-Hartree-Fock 
calculations [32,33] at intermediate Fermi momenta 0.5 fm~ < fcp < 1-5 fm~ . 

Finally, we assume that the interacting particles are on the Fermi surface. 
The momentum variables q, q' and p + p' then form an orthogonal set of basis 
vectors, which follows straightforwardly from the conditions |p ± q/2| = kp 
and |p'±q/2| = kp. Our geometry is shown in Fig. 3, where 9^ and 9^' denote 
the angles between p + q/2, p — q/2 and p + q/2, p' + q/2 respectively and 
are related to the momentum transfer by means off^ 



o = 2 /cp sin — 
^ 2 

9 I 

q' = 2kp sin — 
^ 2 



, /C0s6'q + C0s6'q' 

|p + p I = 2 /cp W ^— ^, 



(12) 
(13) 

(14) 



The orthogonality of the coordinate system is very convenient for the nu- 
merical implementation, in particular for the evaluating the arguments of the 
scattering amplitude in the beta functions. Given this geometry, the scatter- 



^ Our definition of ^q/ differs from the angle between the momenta p and p' 
employed by Ainsworth and Bedell [29]. The relation between the two angles is 
given by 2 cos 6'q' = 1 — cos 0q + (1 + cos 9q) cos ^l- 
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ing amplitude depends only on the magnitude of the two momentum transfers 
q and q' . Under the restriction that the momenta of the quasiparticles are 
constrained to the Fermi surface, we have (f + g'^ ^ 4A;p or equivalently 
COS^q + cosdqi ^ 0. 

The beta functions involve scattering to fast particle-hole intermediate states 
away from the Fermi surface. In the coordinate system shown in Fig. 3 one 
finds 

Q 

P = /cfCOS^P (15) 

p' — kp cos ^ p'. (16) 

Using these relations we obtain, for small q, where p- p' ~ cos6'q', the approx- 
imate expression 

|p + p'|=2A:pcos^cos^, (17) 

which allows us to smoothly extend the RG flow to q^ -\- q'^ > 4/i;|, where 
Eq. (14) is not useful. 



3 The forward scattering limit in the RG and the inclusion of spin 



At the one-loop level, the solution of the RG flow in the particle-hole channels 
is particularly transparent and involves only the effective scattering amplitude. 
The quasiparticle interaction at arbitrary cutoff is simply given by the g = 
values of the scattering amplitude, 

A/(^ = 0,q';A) = Aa(^ = 0,q';A). (18) 



In Fermi liquid theory, the singularity in the ZS channel can be used to elim- 
inate the quasiparticle-quasihole contributions in this channel by taking the 
limit q/uj — > 0. For an energy- independent RG flow, the role of uj is taken 
over by the cutoff. In the RG equation, Eq. (7), the ZS contribution can be 
explicitly eliminated by setting q to zero (see also Fig. 2, where the contri- 
butions from a fast particle and a fast hole are absent for q — iS). Therefore, 
the forward scattering limits in Fermi Liquid Theory (FLT) are in the RG 
approach at zero temperature obtained by 

/FLT(q') = lmia(9 = 0,q';A) (19) 
aFLT(q') = lima(q, q'; A = 0). (20) 
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This is analogous to the q/u ^ and uj/q ^ hmits in Fermi hquid theory. 
The forward scattering amphtude and the Fermi hquid quasiparticle interac- 
tion are obtained as stable solutions of the RG equation, Eq. (7), as the cutoff 
is taken to zero. We note that the beta functions explicitly depend on the cutoff 

and the momentum transfer in the particular channel. In fact, A ^'^(^j ^'i ^) 
is proportional to the cutoff itself, and thus the beta functions vanish by con- 
struction as the cutoff is taken to zero. From the definition of the quasiparticle 
interaction, Eq. (19), we observe that the Fermi liquid parameters depend in 
a nontrivial fashion on the effective scattering amplitude. 

The inclusion of spin is straightforward. The RG equation for the scattering 
amplitude, a(q, q'; A) -|- 6(q, q'; A) cr • tr', is given by 

2A) {/3zs[a,ci,A]+Pzs[b,^,A]cT-a'} 

2A){i/3zs'[«,q',A] + ^/3zs'[&,q',A]} 

2A) /3zs'[«, q', A] - ^ Mb, q'. A]} a ■ a'. 

(21) 

where the coefficients in Eq. (21) are the recoupling coefficients generated by 
the spin exchange operator (see [8]). In [34], the RG flow is illustrated for a 
schematic model of spin polarized fermions. 



A — (a + ha-a')=Q{q- 

- ©(?' - 

- Q{q' - 



4 Results 



We initialize the RG flow with the unique low momentum nucleon-nucleon in- 
teraction Viow k of Bogner et al. [13,14]. This is derived from realistic nucleon- 
nucleon potentials, such as the Bonn, Paris, Argonne and Nijmegen interac- 
tions, by integrating out the high momentum modes while preserving the half- 
on-shell T matrix and consequently the experimental phase shifts. In Viowk 
the short-range repulsion has been tamed. Therefore, it is possible to per- 
form many-body calculations directly with Viow k without first computing the 
Brueckner G matrix. In fact, for Ay^^^^ ~ /cf, ^owk exhibits many of the 
salient features of the G matrix [8] . 

In the case of a static, cutoff independent z^-^ factor, a minimal in medium cor- 
rection to the direct interaction is, as discussed above, included by multiplying 
^lowk by from the start of the RG. After projecting out the non-central 
components, the direct and exchange contribution are in a partial wave basis 
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Fig. 4. The comparison of the RG solution for the forward scattering amphtude 
^FLT (cos 6qi ) and .Bflt (cos 6q' ) (points) to that obtained from the corresponding 
quasiparticle interaction (The sohd hnes include the Fermi liquid parameters up 
to I ^ 5.). The solutions given are for the adaptive z^p factor and a density cor- 
responding to kp = 0.6 fm~^ (dots), fcp = 1.0 fm~^ (squares) and kp = 1.35 fm~^ 
(triangles). Similar agreement is found at other densities. The average momentum 
grid spacing is Aq/kp = Aq'/kp = 0.033. 



given by 



S 



^(q,q';Ao) S) = zl(Ao) N(Ao) 



2S 



X 



|q - q| 



ISJT^l 



-E(2J + 1) (l - (-1)'+^+^) 

q' + q| 



low k 



ISJT^l 



iq + qj ■ (q - qj 
|q' + q| |q' - q| 



(22) 



where we set T = 1 for neutron matter and Pi denotes the Legendre polynomial 
of degree I. We include all partial waves up to J ^ 6 in the driving term. 
Furthermore, we have introduced the dimensionless scattering amplitude A 
by means of 



^(q, q'; A) = A(q, q'; A) + S(q, q'; K) a ■ a' 



(23) 



^ (A) iV(A) ( a(q, q'; A) + 6(q, c(- K) a ■ a'). 



In the scattering geometry discussed in the previous section, q and q' are 
orthogonal, and thus the magnitude of the relative momentum in the initial 
and final states are both given by 



|q^±q| ^ ^34) 
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Fig. 5. The ratio of the direct to the full effective mass m*{Ao)/m*(A = 0), which 
we use as a running z^p factor in the RG equation, versus Fermi momentum kp. 

Let us now discuss the choice of the cutoff in the direct interaction, Ay^^^y.- 
In [8] we argued that for the scattering of two particles with total momentum 
p + p' = relative to the Fermi sea, the choice Av,^^ k ~ accounts for the 
Pauli blocking of the intermediate states. For non-zero total momenta, the 
effect of Pauli blocking is on the average more severe. The Pauli blocking fac- 
tor, averaged over the total momentum of two particles on the Fermi surface, 
corresponds to an angular averaged cutoff of ~ 1.5 kp. We choose the cut- 
off for the low momentum nucleon-nucleon interaction at Avi^^^ = fep-F^ 
This choice has two additional somewhat technical advantages. First, inspec- 
tion of the flow equations, Eqs. (4,5), reveals that momentum transfers up to 
are needed to scatter to intermediate states away from the 
Fermi surface. Since V^ow k is defined only for momenta less than the cutoff, one 
would like to have fcmax < ^Vi^^ k ^ consistent description. Second, a cutoff 
^Mow k ~ kp allows one to define the scattering amplitude continuously on 
the momentum grid q,q' ^2 kp. 

The numerical procedure for solving the RG equations, Eq. (7) or Eq. (21), 
may be checked in the following way. First, the antisymmetry of the effective 
scattering amplitude is by construction preserved under the RG 



Furthermore, the Fermi liquid quasiparticle interaction and the forward scat- 
tering amplitude, Eqs. (19,20), must satisfy general relations of Fermi liquid 
theory. The Fermi liquid parameters are given by the projection of the quasi- 



Since the cutoff dependence of V|ow k is very weak for A > 0.7 fm ^ [14], the exact 
value of Avi^^ is not decisive. 



d 

dA 



{S=l I ^(q, q'; A) + ^(q', q; A) | 5 = 1 ) = 0. 



(25) 
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particle interaction on Legendre polynomials 

^FLTlcOS^qO = zliO) 7V(0)/FLT(q') = E(^' + ^1 ^ ' ^/(cOS^qO, (26) 

where 6^' is the angle between the quasiparticle momenta p and p' in the limit 
q = 0. Similarly, the forward scattering amplitude is expanded in Legendre 
polynomials 

^FLT (cos ^q' ) = E (^^ + ^1 (7 ■ a') Pi (cos 9^^). (27) 

The forward scattering amplitude is related to the Fermi liquid parameters 
through [3] 

A, - 

1 + Fi/{21 + 1) 

P _ Gi 
' l + Gi/i2l + l) 

These relations are fulfilled by our numerical solution with an accuracy of the 
order of Aq/kp and Aq'/kp, where Aq and Aq' are the momentum grid spac- 
ings in q and q' respectively In Fig. 4 we compare the numerical RG solution 
for the forward scattering amplitude to that obtained from the corresponding 
quasiparticle interaction using these relations, Eqs. (28,29). The antisymmetry 
of the scattering amplitude is preserved with machine precision. 

We start the discussion of the numerical results by presenting the change of 
the effective mass, i.e., the ratio of the effective mass obtained from the di- 
rect interaction to the full effective mass. The ratio m*(Ao)/m*(A) is taken 
as an approximation to the running factor in the RG fiow. The final 
value m*(Ao)/m*(A = 0) is shown in Fig. 5. For Fermi momenta 1.0 fm^^ < 
kp < 1.5 fm~^, this approximation agrees well with the values obtained in 
Brueckner-Hartree-Fock calculations including rearrangement terms, as re- 
ported by Zuo et al. [32]. At lower densities, 0.5 fm~^ < < 1.0 fm~^, the 
agreement is qualitatively good, but we find a somewhat smaller value of the 
Zk-p factor compared with the results of Baldo and Grasso [33]. We notice that, 
for densities below k-p = 0.5 fm^^, this procedure yields a quasiparticle strength 
slightly larger than one, while above 1.5 fm""^, the resulting Zk^ increases with 
density. The increase of Zkp with density cannot be excluded from first princi- 
ples, but seems unlikely. Thus, we conclude that our simple procedure seems to 
work reasonably well in the density range 0.5 fm^^ < kp < 1.5 fm^^. However, 
outside this range, corrections due to e.g., rearrangement terms, are probably 
needed. In this paper, our emphasis is on particle-hole polarization effects. 



(28) 
(29) 
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4-1 The Fermi liquid parameters and the scattering amplitude at finite q 

In Fig. 6 we show our results for the / = and / = 1 Fermi hquid parameters 
as functions of the Fermi momentum. Since the RG flow is antisymmetric 
by construction, the Fermi liquid parameters obey the Pauli principle sum 
rules [3,30]. The low momentum interaction drives the flow in the ZS' channel 
for the quasiparticle interaction, and thus our results include effects of the 
induced interaction of Babu and Brown. A further advantage of the RG 
approach is that the calculations can be performed without truncating the 
expansion of the quasiparticle interaction, Eq. (26), at some I. 

Generally, we find a very good agreement between our resultsP^ and the ones 
obtained using the polarization potential model by Ainsworth, Wambach and 
Pines [16,17]. There are minor differences in the value of the effective mass, 
which is treated self-consistently in the RG approach. We note that, in the 
density range 0.6 fm~ < kp < 1.4 fm~ , we find that the effective mass at 
the Fermi surface exceeds unity. The quasiparticle interaction was also calcu- 
lated previously using the induced interaction in [27,35,36]. In these papers, 
the value for the / = spin-dependent parameter Gq is in good agreement 
with ours, while for the spin-independent Fq, there are differences. We stress 
the important role of the large Go for the induced interaction. This Landau 
parameter causes the strong spin-density correlations, which in turn enhance 
the Landau parameter Fq and consequently the incompressibility of neutron 
matter. 

The results can be qualitatively understood by inspecting the RG equation 
for the quasiparticle interaction, Eq. (21) for q = 0, 

Q{q' - 2A) {1 /3zs'[a, q', A] + ^ Mb, q'. A]} (30) 

e(g'-2A) {^/5zs'[a,q',A] - 1 /?zs'[&, q', A]}. (31) 

In this qualitative argument we can neglect Fermi liquid parameters with 1^1. 
At a typical Fermi momentum kp = 1.0 fm~^, we observe that the initial Fq 
and Go are similar in absolute value, |Fo| ~ |Go| ~ 0.8. Consequently, there 



By analyzing the RG equations in [8], one finds that the one-loop RG generates 
a subset of the particle-hole parquet diagrams. The diagrams that are missed cor- 
respond to terms, where the renormalization is due to integrating out internal lines 
of the vertex functions 7 in Fig. 1. In the parquet equations, they correspond to 
certain diagrams involving left or right vertex corrections in the direct or exchange 
particle-hole channel. 

^■^ We note that at large Fermi momenta a smaller z^p factor, as discussed above, 
would lead to an overall decrease of the magnitude of the Fermi liquid parameters. 



A-a(.^0,q';A)^- 
A^% = 0,q';A) = - 
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Fig. 6. The 1 = and I = 1 Fermi liquid parameters versus the Fermi momentum kp. 
The dots denote the direct contribution only (z^p = 1, but including the effective 
mass in the density of states), whereas the squares (constant Zkp factor) and the 
triangles (adaptive Zkp factor) are calculated from the full RG solution. The results 
of Wambach et al. [17] are given for comparison as dashed lines. 

is a cancellation between the contributions due to the spin-independent and 
spin-dependent parameters in Eq. (31), while in Eq. (30) both contributions 
are repulsive. Thus, one expects a relatively small effect of the RG flow on Go 
and a substantial renormalization of Fq, in agreement with our results. 

The approximate relationship for the z^^ factor, Eq. (11), is reasonable, if 
there is a strong frequency dependence of the polarization contributions to the 
nucleon self-energy. This may be the case if one of the Fermi liquid parameters 
is strongly attractive [37]. In the density range where we obtain a reasonable 
Zk-p factor, Fq is indeed most negative. 



The RG approach enables us to compute the scattering amplitude for general 
scattering processes on the Fermi surface, without making further assump- 
tions for the dependence upon the particle-hole momentum transfers q and 
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Fig. 7. The RG solution for the scattering amphtude on the Fermi surface ^(q, q') 
and -B(q, q') versus cos 9qi for different cuts at finite momentum transfer q = 0.5, 1.0 
and V^kp- The RG flow is solved with the adaptive Zk^, factor for different densities, 
= 0.6 fm~^ (dots), fcp = 1.0 fm~^ (squares) and kp = 1.35 fm~^ (triangles). 

q'. In Fig. 7, we show the full scattering amplitude for different values of the 
momentum transfer q and different densities. We recall that for scattering on 
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the Fermi surface the exchange momenta are restricted to + g'^ ^ 4 /c|. The 
scattering amphtude at finite momentum transfer is of great interest for calcu- 
lating transport processes and for assessing the consequences of superfluidity 
in neutron stars. 



4-2 The superfluid gap 



It has been discussed by several authors that polarization effects, which are 
necessary in order to satisfy the Pauli principle in microscopic calculations, 
have a profound influence on superfluidity in neutron stars [16,17,38,39]. A re- 
duction of the pairing gap due to polarization effects was flrst found by Clark 
et al. [38]. Subsequently, Ainsworth, Wambach and Pines explored particle- 
hole polarization effects in detail, using the framework of the induced interac- 
tion [5] combined with the polarization potential model [16,17]. They found 
a large reduction of the ^Sq superfluid gap from a maximum gap of 5 MeV to 
0.9 MeV, when the induced interaction was included. A reduction to similar 
magnitude was obtained by Chen et al. in a variational calculation, when CBF 
correlations were included [40,41]. Recently, Schulze et al. also found a sub- 
stantial reduction of the gap by employing a phase-space averaging method 
to extrapolate the scattering amplitude to back-to-back kinematics [36] . 

Nevertheless, when one compares these calculations in more detail, one finds 
that the maximum gap and the density dependence vary appreciably among 
the different approaches (for a review on superfiuidity of neutron matter and a 
compilation of the results, see the recent notes by Lombardo and Schulze [42]). 
Consequently, further studies are needed in order to improve the microscopic 
calculations and in order to make reliable predictions of superfiuid properties of 
neutron stars. The RG approach offers a systematic framework for computing 
the quasiparticle scattering amplitude, which is needed in the calculation of the 
superfiuid gap. In this paper we have presented a straightforward application 
of the RG method, where both momentum scales, q and q' , are treated on the 
same footing, and consequently the antisymmetry of the scattering amplitude 
is preserved. Moreover, the RG method can provide a unifying and transparent 
framework for the particle-particle and particle-hole effects. Finally, by using 
the unique low momentum interaction Viow k as the direct interaction, much 
of the model dependence is removed. 

We estimate the ^Sq pairing gap using weak coupling BCS theory. The pairing 
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matrix element averaged over the Fermi surface is given byP^ 



(A) = ^ £ dcos9^{A{ q, ^Akl - )-3B{q, ^ Ak^ - q^ ) } , (32) 

where the momenta for back-to-back scattering are constrained by q^ + q'"^ = 
4kp. An estimate of the superfluid gap is then given by the weak couphng 
formula [44] 

A = 2eF exp(^^^, (33) 

where ep = kp/2m* denotes the Fermi energy and we use the same prefactor 
as Ainsworth et al. [16,17]. For liquid ^He, the weak coupling BCS results for 
the critical temperature are in remarkable agreement with experiment when 
experimental Fermi liquid parameters are used as input [44,45] . For the direct 
interaction, the pairing gap, Eq. (33), with a free single particle spectrum is 
given by 

A = ^expf ^ --Y (34) 

m '\2kFmVio^y,{kF,kr,A = V2kr,^So)J 



In Fig. 8 we show the superfluid gap computed with the direct term only, 
Eq. (22), as well as that obtained with the RG solution for the scattering 
amplitude. To facilitate the comparison with other calculations, the gap ob- 
tained with the direct interaction is given for a free single particle spectrum, 
i.e., m*/m = Zkp = 1. The full solution of the BCS equation using various 
free nucleon-nucleon interactions yields a maximum gap of about 3 MeV at a 
Fermi momentum of 0.8 — 0.9 fm~^ (for an overview of the results see Fig. 2 
in [42]; see also [46,47,48]). The agreement of the weak coupling result us- 
ing Viowk, Eq. (33), with these calculations is remarkably good. The reason 
for this success is presumably that in V^ow k, the strong short-range repulsive 
core is eliminated, whereas the phase shifts are preserved. For such a "weak" 
interaction, the weak coupling treatment is expected to work. 

In the full RG calculation, including particle-hole polarization effects, we flnd 
a much smaller gap; the maximum gap is 0.8 MeV at a Fermi momentum of 
kp = 0.8 fm~^. Here, both the effective mass and the quasiparticle strength 
zl^ are taken into account in the quasiparticle scattering amplitude. The gap 
is in good agreement with that obtained by Wambach et al. [17], which is 

-"^'^ The factor 1/8 in Eq. (32) consists of a factor 1/2 for the average over cos^q, 
another factor 1/2 for the reduction of the full density of states to that of one 
spin species. The remaining factor 1/2 accounts for the normalization of the anti- 
symmetrized contact interaction employed in the standard derivation of the weak 
couphng result (see e.g., [43]). 
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Fig. 8. The Sq superfluid gap versus the Fermi momentum kp. The dots denote 
the direct contribution only {m* /m = z^p = 1), whereas the squares (constant 
Zfcp factor) and the triangles (adaptive factor) are calculated from the full RG 
solution. In comparison we show the values of the gap calculated by solving the 
full BCS equation (dashed line), taken from [42]. The superfluid gap including par- 
ticle-hole polarization effects obtained by Wambach et al. [17] is shown as a solid 
line. 

shown for comparison in Fig. 8. The results obtained in other calculations, 
which include polarization effects, are compiled in Fig. 7 of [42]. 

We observe that the full gap vanishes at k-p = 1.45 fm~^ both for a constant 
and an adaptive z^p factor. Moreover, at the lowest Fermi momentum we 
considered, kp = 0.2 fm~ , the ratio of the full to the direct gap is given 
by A/Aq = 0.43, which is close to the universal low density limit A/Aq = 
(4e)-V3 ~ 0.45 [49].[3 

In the direct interaction, Viowk, particle-particle ladder diagrams are included 
for large relative momenta, k > Ayj^^j^. In fact, we argued that the cutoff 
accounts for the Pauli blocking of intermediate particle states. This raises the 
question of double counting, since the BCS equation also sums particle-particle 
ladders. However, since the gap is generated by low-lying two-particle and two- 
hole states, which are excluded in Viow k due to the cutoff at V^kp, the double 



We note that the low density limit for the ratio of the full to the direct S-wave 
gap is independent of the scattering length. While the strict low density expansion 
in kpas is not relevant to neutron star matter, since its validity is restricted to 
densities below neutron drip, an expansion in an effective scattering length kpa^s 
may be valid up to much higher densities. As argued in [8], Vjow k(0, 0; Ay,^^ ^ ~ kp) 
can be interpreted as the effective scattering length in the particle-particle channel. 
In the isotriplet channel V|owk(0)0) ~ 2.0 fm in the cutoff independent region. If 
one adopts the effective scattering length, it is not surprising to find a ratio close 
to the universal one at kp = 0.2 fm~^, since kp <C l/jaeffj ~ 0.5 fm~^. 
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counting problem is minimal. Furthermore, because the cutoff dependence of 
^ow k is very weak in the density range relevant for superfluidity, any residual 
double counting, due e.g., to states above the cutoff that contribute to the 
solution of the BCS gap equation, would not affect our results significantly. 

We note that by using a density dependent cutoff, we recover the correct 
low density limit, where the effective neutron- neutron interaction is given by 
the free T matrix. This follows because the low momentum interaction Viow k 
approaches the T matrix for momenta below the inverse S-wave scattering 
length, i.e., A^^^j^ ~ kp <^ l/|as|, as shown in Fig. 5 of [14]. In particular, 
Viow k{kF,kp) as/m, as kp 0,|2^ and we recover the well known low den- 
sity results in Eq. (34), where the gap is given by the phase shifts [50]. However, 
due to the large neutron-neutron scattering length, the low density expansion 
is applicable only for Fermi momenta below kp -C l/|as| = 0.05 fm~^, which 
corresponds to extremely small densities, well below neutron drip. Since neu- 
tron star matter at such densities is in a solid phase, the strict low density 
expansion in kpa^ is not relevant for the physics of neutron stars. 

Finally, we remark that the physical pairing gap is given by Zkp A. Thus, our 
results should be multiplied by the quasiparticle strength given in Fig. 5.p^ 
We stress that the physical effect responsible for the strong reduction of the 
gap is the polarization by the many-body medium. A detailed analysis of 
the influence of self-energy contributions on the superfluid gap can be found 
in [33,51,52]. 

Particle-hole polarization may have a strong effect also on the color supercon- 
ducting gap in high density QCD. We plan to explore these effects using the 
methods presented here in future studies. 

4-3 The equation of state of neutron matter 

In contrast to nuclear matter, neutron matter does not saturate. It is stable 
only if the positive pressure is balanced by an external force, as e.g., gravitation 
in a neutron star. Therefore, it is possible to compute the equation of state 
of neutron matter using the Fermi liquid parameters presented in Section 4.1. 
This method, which accounts for polarization effects on the equation of state, 
was first employed by Kallman [15]. 

The speed of sound cs of a Fermi liquid is given by the incompressibility K [53] 
In [14] conventional units are used where hi^/m = 1. 

Given the uncertainty in the pre-exponential factor of Eq. (33) , this may seem 
superfluous. However, for a given pref actor, it leads to an additional suppression of 
the gap. 
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^ dp~9m^'' ^^^^ 
where p denotes the mass density and K is the incompressibihty. The latter 
is in Fermi hquid theory given by 

3fci(l + Fo) 

In the density range of interest, the neutrons may to a good approximation 
be treated non-relativistically. Thus, the mass density is given hj p = mn, 
where n = /cp/37r^ is the number density. Now, we can solve the differential 
equation 

dn 9 ^ ^ 

for the equation of state P{n). The incompressibihty is computed at each 
density using the Fermi liquid parameters obtained in the RG approach. The 
energy per neutron is then computed by integrating the thermodynamic rela- 
tion 

P = „^^. (38) 

dn 

This approach requires two initial conditions for the differential equations, 
Eqs. (37,38). One can match our results with a conventional equation of state 
Pc{n) at a certain matching density tt-m- For n > nyi, the resulting equation 
of state P{n) reads 

P(n) = Pin) - Pinu) + PcI^m), (39) 

with a corresponding equation for the energy per particle. We match the pres- 
sure and the energy to the V^ow k Hartree-Fock equation of state at the density 
n = 0.02 fm~^. This is the lowest density given by Akmal et al. [54] We 
stress that the equation of state depends on the matching condition. The en- 
ergy is particularly sensitive, since it requires two constants of integration. 
Below we give a parametrization of the pressure, with the boundary condition 
P{n = 0) = 0, 

P(„) = 79.09 1.000 - 1.177 (j^)"' 

With this parametrization, the reader can match our equation of state to a 
conventional one at low densities. 



At very low densities the Hartree-Fock calculation cannot be trusted due to the 
large neutron-neutron scattering length, as discussed above. 
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Fig. 9. The pressure P and the energy per neutron E/A for pure neutron matter as 
functions of density n. Results are given for the full Fermi liquid parameters (solid 
line) and the direct ones (dashed line) matched to the Vjow k Hartree-Fock equa- 
tion of state at low density, n = 0.02 fm~^. Also given is the Vjow k Hartree-Fock 
contribution to the energy [55] and the pressure calculated by means of Eq. (38) 
(triangles). In comparison, we show the Brueckner-Hartree-Fock equation of state 
of Bao et al. [56] (dots) computed with the non-relativistic Bonn potential and the 
equation of state of Akmal et al. [54] (squares) from a variational calculation using 
the Argonne vis interaction supplemented by three-body forces and the boost cor- 
rections. For a comparison with other equations of state, see also Fig. 1 in Lattimer 
and Prakash [57]. 

The resulting pressure and the energy per neutron are shown in Fig. 9. We note 
that the Hartree-Fock energy obtained with Viovj k and the energy obtained 
from the direct Landau parameters do not coincide. This is solely due to 
the density dependence of the cutoff A^^^^, which gives rise to additional 
contributions to the effective interaction. However, since a major part of these 
contributions correspond to particle- hole propagation in the ZS' channel, they 
are to a large extent part of the induced interaction and should not be included 
in the direct term. Consequently, the difference between the Hartree-Fock 
energy and that obtained using the full induced interaction is due to higher 
order particle-hole polarization effects [10]. 

Relative to the l^iow k Hartree-Fock calculation [55] there is an increase in 
pressure and energy below nuclear matter density. Compared to other non- 
relativistic equations of state [56,54], the pressure and energy is somewhat 
larger at low densities At densities above 0.11 fm~^, corresponding to Fermi 
momenta kp > 1.5 fm~^, the resulting z^p factor is larger than that found 
in Brueckner-Hartree-Fock calculations, as noted above. A smaller z^p factor 
would reduce the attractive contribution from Fq and thus lead to an increase 



° The equation of state of Akmal et al. is significantly stiffer at higher densities due 
to three-body forces. 
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of the incompressibility and consequently also of the pressure and the energy 
for n > 0.11 fm^^ in our calculation. We conclude that polarization effects 
possibly have an important impact on the equation of state, which should be 
included in a consistent treatment of neutron star superfluidity and neutron 
star structure. 



5 Summciry and conclusions 

In this work, we have developed a practicable framework for renormalization 
group calculations of strongly interacting Fermi systems. We employed the 
RG approach to compute the quasiparticle interaction and the scattering am- 
plitude on the Fermi surface for neutron matter, which for many purposes is 
a good approximation to neutron star matter, where there is a small admix- 
ture of protons. Our approach follows the ideas of Shankar in applying RG 
techniques to the interacting fermion problem [6,19]. 

RG methods used to solve multi-channel problems have several advantages. 
First, all channels and the corresponding momentum scales can be treated on 
equal footing in a systematic way, which can accommodate fundamental sym- 
metries, such as the Pauli principle. In this study, we considered the RG flow 
in the particle-hole channels at the one-loop level. An expansion in momen- 
tum transfers was performed, which enabled us to compute the beta functions 
analytically. We note that it is fairly straightforward to include the particle- 
particle channel in the RG equations for the scattering amplitude and the 
quasiparticle interaction. This opens the possibility to explore the interfer- 
ence of the particle-hole and particle-particle channels. A calculation of this 
kind corresponds to solving the "parquet equations at the one-loop order" . The 
expansion in momentum transfers can be removed by computing the complete 
particle- hole phase space [12]. 

In contrast to the traditional approach to the many-body problem, the RG 
adapts the starting vacuum interaction to include the in medium effects due to 
the fastest particle-hole excitations. In the subsequent iteration, this renormal- 
ized interaction is modifled to account for the next fastest modes. Therefore, 
at every step, the flow of the effective interaction is expanded around the cur- 
rent system, which highlights the efficacy of the approach. Moreover, the RG 
flow is physically very transparent, since the effects of the high momentum 
modes are readily tractable. We believe that the RG method is a promising 
tool for studying a wide range of nuclear many-body problems. 

At the one-loop level and at zero temperature the RG equations are partic- 
ularly transparent. We found an ambiguous long- wavelength limit, analogous 
to that appearing in the original treatment of Landau, which can be used to 



25 



distinguish between the quasiparticle interaction and the forward scattering 
amphtudc. This imphes that, at the one-loop level, it is necessary to solve 
only the RG equation for the scattering amplitude. The quasiparticle interac- 
tion is then obtained in the corresponding hmit. A test for the consistency of 
the solution is provided by general relations between the forward scattering 
amplitude and the Fermi liquid parameters. 

We solved the RG equations in the particle-hole channels for neutron matter, 
where tensor interactions are absent in relative S-states. The RG flow includes 
the exchange channel, and thus it builds up the polarization effects of the 
induced interaction of Babu and Brown. Our results for neutron matter arc 
very encouraging and provide strong motivation for further developments of 
the RG approach to the nuclear many-body problem as well as for condensed 
matter systems. 

In the calculation, polarization and self-energy effects have been treated self- 
consistently. The latter have been included in a simple approximation, which 
nevertheless led to good quantitative results at intermediate densities. 

Results for the Fermi liquid parameters and the scattering amplitude for the 
density range of interest to neutron stars were presented. The Fermi liquid 
parameters satisfy the Pauli principle sum rules by construction. Moreover, 
polarization effects lead to an enhancement in the pressure and the energy 
per neutron below nuclear matter density. As a first application, which probes 
the momentum dependence of the scattering amplitude, we computed the ^Sq 
superfluid gap using weak coupling BCS theory. A reduction from a direct 
gap of 3.3 MeV to 0.8 MeV was found. Our results confirm the importance 
of particle-hole polarization effects for superfluidity. A further application of 
the scattering amplitude at finite momentum transfers is in calculations of 
transport processes. 

For a similar analysis of nuclear matter, it is necessary to extend the flow 
equations to incorporate tensor interactions. These are important for a com- 
plete discussion of nuclear matter, and we expect rather large renormalization 
effects from the arguments given in [8] . Moreover, tensor correlations and the 
tensor Fermi liquid parameters play an important role in the physics of dense 
matter [58] . As noted in the introduction, tensor forces are less pronounced in 
neutron matter, since they are not operative in relative S-states. However, for 
speciflc quantities the tensor force may be of crucial importance. 

To date, polarization effects on the ^P2-^F2 pairing of neutrons have not been 

addressed and a quantitative analysis is very much needed. Pethick and Raven- 
hall argue that both density and spin-density fluctuations would increase the 
gap in this state and thus possibly extend the superfluid regime to lower densi- 
ties [39] . Once the RG approach is adapted to allow for tensor forces, it should 
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be straightforward to address this problem [59]. 
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A Derivation of the Fermi liquid theory results from the ZS chan- 
nel beta function 



For a one channel problem, such as the interacting Fermi system restricted to 
the ZS channel, the one-loop RG is equivalent to solving the Bethe-Salpeter 
equation in that channel. We define the fast particle-hole Lindhard function 
at a; = by integrating out the fast modes 

x{q,A)^g / 7^-^^^^^ / • (A.l) 

For q ^ kp, we expand and write the particle-hole phase space explicitly at 
zero temperature 



x{q,A) = 9 j 



X 



(27r)^ v^qcosO 

e(A;F - A - (/ + gcos^/2)) Q{{p" - gcos^/2) - (fcp + A)) 



- e((p" + q cos e/2) - {kp + A)) e(A;F - A - (/ - q cos 9/2)) 

(A.2) 

where the occupation numbers are automatically taken care of by the func- 
tions that define the border between the fast and slow modes and 9 denotes 
the angle between p" and q. We rewrite the product of theta functions using 
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X 



e(Ai - x) e{x - A2) = (e(x - A2) - e{x - Ai)) e(Ai - A2). This leads to 

d^p" 1 

(27r)^ vpgcosO 

{ e(p" - g COS e/2 - A;f - A) - e(p" + g COS ^/2 - A;f + A) } 
0(-g cos ^ - 2 A) - { &{p" + q COS ^/2 - (A;f + A)) 



X 



X 



e(p" -qcose/2-kF + A)} e(q cos e - 2A) 



(A.3) 



When we differentiate with respect to the cutoff, the contributions from the 
multiphcative theta functions 0(±g cos ^— 2A) vanish, since for q cos 6/2 = ±A 
the two remaining theta functions cancel. Moreover, q <^ kp implies A <^ kp, 
which leads to 



dA 



1 



(27r)^ vpqcosO 
e{-2A-qcos9) 

2 gm*k-p 



2 5{p"-kp) 



e(gcos^-2A) 



(A.4) 



q 2 7r^ J2A/q<|cos6»|<l An 



dClr,,, signfcos^j 

P ^ ^ e(g-2A). (A.5) 



cos^ 



This leads to the RG equation, Eq. (7). Under restriction to the ZS channel, 
we find after projecting the scattering amplitude on Legendre polynomials 



a = ^ai(A)P,(cos^q/) 



I 



daijA) 
dA 



27r2 



q I 21 + 1 



(A.6) 
(A.7) 



where we have set Zkp = 1 for simplicity. We propose to think of q as an 
external scale, since it must be kept finite and may only be set to zero after 
integrating the RG equation. The solution of the RG equation in one channel, 
Eq. (A.7), is obtained by integrating from the initial condition a/(A = q/2) — 



,(0) 



to A = 0. We find 



,(0) 



ai 



(0) 



1 + 



g m^kp tti 
27r2 21 + 1 



(A.8) 



which are the standard Fermi liquid theory relations among the quasiparti- 
cle interaction and the forward scattering amplitude after identifying Ai — 



gm*k^ 



ai and Fi 



af'^ . The RG initial condition being the Fermi liq- 



2^2 - ' 27r2 

uid parameters is understood, since we assume that the flow in the ZS' channel 
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has been carried out to build up the quasiparticle interaction or equivalently 
the Fermi hquid parameters are taken from experimental measurements. 



References 

[1] L.D. Landau, Sov. Phys. JETP 3 (1957) 920. 
[2] L.D. Landau, Sov. Phys. JETP 5 (1957) 101. 
[3] L.D. Landau, Sov. Phys. JETP 8 (1959) 70. 

[4] G. Baym and C.J. Pethick, Landau Fermi Liquid Theory: Concepts and 
Applications, Wiley, New York, 1991. 

[5] S. Babu and G.E. Brown, Ann. Phys. 78 (1973) 1. 

[6] R. Shankar, Rev. Mod. Phys. 66 (1994) 129. 

[7] J. Polchinski, in Proceedings of the 1992 Theoretical Advanced Studies Institute 
in Elementary Particle Physics, Eds. J. Harvey and J. Polchinski, World 
Scientific, Singapore, 1993, |hep-th/9210046 . 



[8] A. Schwenk, G.E. Brown and B. Friman, Nucl. Phys. A703 (2002) 745. 
[9] K.S. Bedeh and T.L. Ainsworth, Phys. Lett. 102A (1984) 49. 



[lo; 

[11 
[12 

[13; 

[14 



[15 

[16; 
[ir 

[is; 

[19 

[2o; 

[21 
[22; 



O. Sjoberg, Ann. Phys. 78 (1973) 39. 

S.-O. Backman, G.E. Brown and J.A. Niskanen, Phys. Rept. 124 (1985) 1. 
V. Brindejonc and B. Friman, in preparation. 

S.K. Bogner, T.T.S. Kuo and L. Coraggio, Nucl. Phys. A684 (2001) 432c. 



S.K. Bogner, T.T.S. Kuo, A. Schwenk, D.R. Entem and R. Machleidt, |nucl 



th/0108041 



C.-G. Kallman, Nucl. Phys. A213 (1973) 461. 

T.L. Ainsworth, J. Wambach and D. Pines, Phys. Lett. B222 (1989) 173. 

J. Wambach, T.L. Ainsworth and D. Pines, Nucl. Phys. A555 (1993) 128. 

K.G. Wilson and J.B. Kogut, Phys. Rept. 12 (1974) 75. 

R. Shankar, Physica A177 (1991) 530. 

N. Dupuis, Eur. Phys. J. B3 (1998) 315. 

G. Benfatto and G. GaUavotti, Phys. Rev. B42 (1990) 9967. 

G. Benfatto and G. GaUavotti, J. Stat. Phys. 59 (1990) 541. 



29 



[23] J. Feldman, D. Lehmann, H. Knorrer and E. Trubowitz, in Constructive 
Physics, Springer Lectures Notes in Physics, Ed. V. Rivasseau, Vol. 446 (1995) 
267. 

[24] N. Dupuis and G.Y. Chitov, Phys. Rev. B54 (1996) 3040. 
[25] A.C. Hewson, Adv. Phys. 43 (1994) 543. 

[26] S.-O. Backman, O. Sjoberg and A.D. Jackson, Nucl. Phys. A321 (1979) 10. 

[27] S.-O. Backman, C.-G. Kallman and O. Sjoberg, Phys. Lett. B43 (1973) 263. 

[28] T.L. Ainsworth, K.S. Bedell, G.E. Brown and K.F. Quader, J. Low Temp. Phys. 
50 (1983) 317. 

[29] T.L. Ainsworth and K.S. Bedell, Phys. Rev. B35 (1987) 8425. 

[30] B.L. Friman and A.K. Dhar, Phys. Lett. B85 (1979) 1. 

[31] J.P. Jeukenne, A. Lejeune and C. Mahaux, Phys. Rept. 25 (1976) 83. 

[32] W. Zuo, G. Giansiracusa, U. Lombardo, N. Sandulescu and H.-J. Schulze, Phys. 
Lett. B421 (1998) 1. 

[33] M. Baldo and A. Grasso, Phys. Lett. B485 (2000) 115. 

[34] A. Schwenk, B. Friman, S.K. Bogner, G.E. Brown and T.T.S. Kuo, Talk given 
at 5th Int. Conference Renormalization Group 2002, Tatranska Strba, Slovakia, 
March, 2002, Acta Phys. Slov. 52 (2002) 207. 

[35] A.D. Jackson, E. Krotschek, D.E. Meltzer and R.A. Smith, Nucl. Phys. A386 
(1982) 125. 

[36] H.-J. Schulze, J. Cugnon, A. Lejeune, M. Baldo and U. Lombardo, Phys. Lett. 
B375 (1996) 1. 

[37] K. Bedell, private communication. 

[38] J.W. Clark, C.-G. Kallman, C.-H. Yang and D.A. Chakkalakal, Phys. Lett. B61 
(1976) 331. 

[39] C.J. Pethick and D.G. Ravenhah, Ann. N.Y. Acad. Sci. 647 (1991) 503. 

[40] J.M.C. Chen, J.W. Clark, E. Krotschek and R.A. Smith, Nucl. Phys. A451 
(1986) 509. 

[41] J.M.C. Chen, J.W. Clark, R.D. Dave and V.V. Khodel, Nucl. Phys. A555 
(1993) 59. 

[42] U. Lombardo and H.-J. Schulze, in Physics of Neutron Star Interiors, Springer 
Lecture Notes in Physics, Eds. D. Blaschke, N.K. Glendenning and A. 
Sedrakian, Vol. 578 (2001) 30, |astro-ph/0012209| . 

[43] A. A. Abrikosov, L.P. Gor'kov and I.E. Dzyaloshinski, Methods of Quantum 
Field Theory in Statistical Physics, Dover, New York, 1963. 



30 



[44 
[45; 

[46; 

[47; 
[48; 
[49; 
[5o; 

[51 
[52 
[53; 

[54; 

[55 
[56; 

[57; 

[58; 



B.R. Patton and A. Zaringhalam, Phys. Lett. A55 (1975) 95. 
K. Bedell and D. Pines, Phys. Lett. A95 (1980) 281. 

M. Baldo, J. Cugnon, A. Lejeune and U. Lombardo, Nucl. Phys. A515 (1990) 
409. 

V.A. Khodel, V.V. Khodel and J.W. Clark, Nucl. Phys. A598 (1996) 390. 

0. Elgar0y and M. Hjorth-Jensen, Phys. Rev. C57 (1998) 1174. 

L.P. Gorkov and T.K. Melik-Barkhudarov, Sov. Phys. JETP 13 (1961) 1018. 

V.J. Emery and A.M. Sessler, Phys. Rev. 119 (1960) 248. 

P. Bozek, Nucl. Phys. A657 (1999) 187. 

U. Lombardo, P. Schuck and W. Zuo, Phys. Rev. C64 (2001) 021301(R). 

P. Nozieres, Theory of Interacting Fermi Systems, Addison- Wesley, Reading, 
Massachusetts, 1997. 

A. Akmal, V.R. Pandharipande and D.G. Ravenhall, Phys. Rev. C58 (1998) 
1804. 

A. Schwenk, Ph.D. Thesis, SUNY Stony Brook, 2002. 

G. Bao, L. Engvik, M. Hjorth-Jensen, E. Osnes and E. 0stgaard, Nucl. Phys. 
A575 (1994) 707. 

J.M. Lattimer and M. Prakash, Astrophys. J. 550 (2001) 426. 

E. Olsson and C.J. Pethick, Talk given at Compact Stars in the QCD Phase 



Diagram, Copenhagen, 2001, |astro-ph/02030Tl| 



[59] A. Schwenk, B. Friman and G.E. Brown, in preparation. 



31 



